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ABSTRACT 


This study is to assess the feasibility of using a laser's 
radiation to remove masses on the order of milligrams from 
a target surface and to form a vapor jet which could then 
serve as the conducting channel in a giant power pulse dis- 
charge. 

A model based on an energy balance is presented which 
predicts the mass removed and the pressure, density and 
temperature of the material as it propagates from the target 
surface as functions of the pulse time and the absorbed power 
density.  Combinations of laser pulse time and power density 
are considered for which the surface temperature remains low 
enough in order to avoid cut-off by a plasma plume in front 
of the target surface. 

A model of adiabatic jet expansion is presented to predict 
the radial expansion of the vapor plume as it propagates in 
the axial direction. Temperature, pressure and density are 
predicted after radial expansion using the model of adiabatic 
jet expansion. 

These models predict, theoretically, that the desired 
masses and limits on vapor plume expansion can be achieved 
using relatively low power densities for the proposed experi- 


ment. 
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I. INTRODUCTION 


Requirements for material and equipment tests at very 
high x-ray radiation fluxes with particular spectral charac- 
teristics have led to the development of giant pulse power 
devices, Johnson [1]. 

The technique of using "exploding wires" or wire arrays 
dates back to 1773 [2]. In this technique the energy from 
a large storage capacitor is discharged first through a pulse 
forming network then through thin metallic wires which pinch 
due to the development of large magnetic fields and eventually 
develop into a high temperature plasma which becomes a source 
of high intensity x-rays. In recent time the wire array has 
been replaced by a gas-puff from a supersonic nozzle. After 
preionization with microwave radiation the discharge through 
Meee gas jet leads to a high temperature pinch discharge. 

The x-ray intensity and spectral distribution depends 
upon the high voltage pulse form, plasma temperature and 
dynamics, mass of the wire array or gas in the gas-puff, 
z-number of the material and the overall geometry of the 
experiment.  Elaborate computer codes, Colombant [3], have 
been written to model the x-ray emission as a function of these 
parameters. In order to obtain the desired spectra and inten- 
sities, these computer codes predict masses on the order of 
milligrams and high z-number materials are required in the 


discharge channel. 





Wire arrays are not available in small enough diameter 
sizes for high z-number materials. As for the gas-puff 
device, materials again are severely limited to gaseous ele- 
ments such as Ar, Kr, and Xe. Because of the shortcomings 
of wire arrays to satisfy the mass and z-number requirements, 
alternate methods of establishing high energy discharge 
channels are of interest. 

The use of high power laser pulses directed at a target 
material has been proposed by Schwirzke, described by Johnson 
[1]. In such a scheme, the target material in the form of a 
solid disc can be mounted on the anode (or cathode) of a high 
voltage pulsed power machine. A laser beam is directed toward 
the target and a certain mass of the target is vaporized, 
figure l. This vapor plume will expand through the evacuated 
chamber toward the cathode. Electrical energy would then 
couple to the vapor plume producing a pinched discharge 
channel with an intense x-ray emission. Preionization of the 
vapor could be done by microwaves as in the gas-puff experi- 
ment or conceivably by the laser itself. The advantages of 
this concept over wire arrays and gas-puffs are great in 
principle as targets could conceivably be made of almost any 
material. 

By controlling the energy, pulse width and focal area of 
the laser along with the choice of material, the x-ray emission 
spectrum could be varied. Combinations of different materials 


could also be used. 
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There are two distinct mechanisms by which a laser can 
remove material from a target. First, the laser energy is 
absorbed by the target to increase the temperature to the 
vaporization temperature. At high power densities this takes 
a very short time, on the order of a few nanoseconds. Since 
the power density is large, a high percentage of ionization 
in the vapor cloud or plume will be reached after a very short 
time. This high ionization level causes absorption of the 
laser energy in the plume thus shielding the target surface 
from the continuing incoming radiation. However, further 
vaporization of the surface via electronic heat conduction 
continues. This mechanism has been investigated by Brooks 
[4]. He used a Q-switched neodymium glass laser with a 5 
Joule output of pulse duration 25 nsec directed at an aluminum 
target with a focal diameter of 0.15 cm. This results ina 
power density of about o>” πο . Measurements showed 
mass removal to be about 0.7 mg. While the mass yield is of 
the desired order, an important feature of this mechanism makes 
it less attractive for the purpose here. The main plasma, 
carrying the bulk of the ionized mass, propagated at mU 
cm/sec. If anode-cathode separation were assumed to be 4 cm, 
the main plasma would reach the anode in 364 nsec. The target 
temperature is so high that vaporization continues for almost 
l0 usec after the laser pulse is stopped. As a result, the 
plasma particles in the plume will expand, thus filling the 


entire cavity before sufficient mass has been removed from the 


target. As the plasma now covers a large volume, energy 





coupling from the high voltage source will become difficult. 
Contamination of and arcing in the gap of the coaxial delay 
line forming the cathode cone will reduce efficiency. 

An alternate approach may be used to overcome this prob- 
lem. If the laser power density is lowered the ionization 
level can be reduced to a point where no significant absorp- 
tion of the laser's radiation takes place in the vapor plume. 
This plume is now a true vapor plume and not a plasma. The 
laser energy is not shielded from the target surface. As the 
surface temperature is now relatively low, the plume propagates 
much slower. The timing of the application of the high voltage 
can be made to correspond to the travel time of the vapor 
plume. The question is whether significant masses, on the 
order of milligrams, can be removed from the target. This is 
the approach that will be studied in this paper. 

There are many parameters controlling the mass removal 
process. The energy delivered by the laser to the material 
surface divides into the reflected portion and the absorbed 
portion. Classically the absorption is dependent on the laser 
radiation wavelength and the conductivity of the target mater- 
ial. In section II , the relevant absorption processes are 
described for various materials of interest and a table of 
absorption is compiled which is used as the basis for the 
absorption calculated. 

The absorbed energy at the surface raises the temperature 
of the material leading eventually to melting and vaporization 


beginning at the surface. The initially steep temperature 
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gradient near the surface leads to rapid heat conduction 
into the material. This process is governed by the heat 
conduction equation. Because of the three phases present 
this is a very complex problem. There are many models of 
this process using various approximations in the literature, 
often employing elaborate computer Programs, In order to 
stay with the purpose of this study, to assess the general 
feasibility of the proposed mechanism, a very simple model of 
the heat conduction and evaporation has been adopted which 
will be discussed in Section III. Its relation to other models 
in the literature are also discussed. 

As the material evaporates from the surface it expands 
both axially and radially. Section IV deals with the 
radial expansion which is described as an adiabatic jet expan- 
sion. Density, temperature and pressure changes in the expan- 


sion process are also described. 
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II. THE ABSORPTION PROCESS 


Materials can be divided into three basic classes, 
conductors, semi-conductors and dielectrics. Dielectric 
materials are not of interest here because they do not con- 
duct well. In order to vaporize these materials, high power 
densities must be used which would result in plasma production 
at the target surface. This is a condition we are trying 
to avoid as these high temperature plasmas lead to surface 
shielding. 

Semi-conducting materials in general are also poor conduc- 
tors. However, at certain wavelengths, there are resonances 
where conduction is quite high. In these regions their thermal 
properties resemble those of metals. Therefore, our considera- 
tions will be confined to the metallic absorption process. 

The interaction of materials with electromagnetic waves 
is formulated on the basis of Maxwell's equations. This 
development is well known and documented [5, 6, 7]. Derivation 
of optical quantities is presented in Appendix A. 

Plots can be made of absorbtivity, À vs. A, the wavelength, 
for several materials using the following expression: 


πετ... (eq. 2-1) 


(ntl) + k’ 


where n and k are the real and imaginary parts of the complex 


EUudex Of refraction. The n and k are both functions of 


p» 





frequency, o, the conductivity, o, and the parameter m*/N 
where m* is the effective mass and N is the free electron 
density of the metal. In general it is very difficult to 
obtain an appropriate value for m*. Much controversy has 
arisen on how to measure effective mass, Seitz [8] lists some 
approximate values which are shown in Table 1l. Harrison and 
Webb [9] review analytical and experimental work which illus- 
trates the difficulty in obtaining appropriate values. This 
problem of choosing a value for m* need not present an obsta- 
cle, however.  Wieting and Schriempf [7] show proof that the 
sensitivity of reflectivity to the parameter ß/N, where 

β Ξ m*/m, y is small. Here m. is the free electron mass. 
Therefore, the selection of values for m* in [8] are justified 
for our purpose. 

Figures 2-9 are plots of n, k, and Á as functions of 
decreasing wavelength using the data from Table I. It is 
interesting that the values for absorbtivity are very low for 
the region where w is ae than the plasma frequency. The 
absorbtivity then increases very sharply at the plasma fre- 
quency for the electrons inthe metal and quickly reaches 100 
percent. The degree of absorption remains relatively constant 
over the infrared region and is small, generally 1-2 percent. 

The question now arises as to whether the calculated 
values for absorbtivity are reasonable. Bennett, Silver and 
Ashley [10] have shown surprisingly close correlations between 


experiment and these classically calculated absorbtivities in 
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ultra high vacuums of ΠΝ. 


torr for aluminum. However, at 
lower values of vacuum the theory breaks down and the absorb- 
tivity is greater. This is significant because in this experi- 
ment the base pressure in the vacuum chamber is Ee τους, 
Bonch-Bruevich [13] also shows experimental results illus- 
trating the change in reflectivity as a function of pulse time 
in copper and steel, figure 10. It must be pointed out, how- 
ever, that these results pertain to conditions different from 
those of interest here. This particular experiment was carried 


out using a laser with a power density on the order of Lows 


Watts/cm^. At this high power level there may be major plasma 
production and the marked decrease in reflectivity at the tar- 
get surface is probably due more to the plume absorption rather 
than the temperature dependent changing of optical properties 
of the target material. The first plateau at 2 usec is probably 
due to the onset of vaporization. One would not expect a much 
lower power level laser to produce the degree of change in 
absorbtivity that Bonch-Bruevich observed. 

In calculating the absorbtivity, dc conductivity was used. 
Conductivity is inversely proportional to the resistivity. 
The resistivity of metals increases nearly linearly with 
temperature in the region from room temperature to melting. 
Figure ll is a plot of absorbtivity as a function of temperature 
for four metals, aluminum, silver, cadmium and zinc, all at 
A = 10.6 um. The underlying assumption in this plot is that 


the liquid and solid phase conductivities are the same. The 
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plot shows increases in absorbtivity on the order of 1-3 
percent. 
Thomas and Musal [14] have proposed using an effective 


value for absorbtivity, Á using the following expression: 


τας. 


s - -lu. -1/2 
A zo nd A 


eff (eq. 2-2) 


where m is the absorbtivity at vaporization temperature and 
A is at room temperature.  Harrach [15] points out that these 
effective values of absorbtivity may be 5-10 times those at 
room temperature and cites the Bonch-Bruevich report, [13], 

as experimental evidence. Again, this could be misleading 

at low power densities. 

Surface preparation plays a large role in the absorbtivity 
of any practical metal surface and is an empirical matter. 
Schriempf [16] states some values for Al, Au, Cu, Ag, Ti, 

Mg and 304-steel. Sandblasting the surface appears to be the 
best technique resulting in up to ten times the abosrbtivi- 
ties of ideal surfaces. He does not state the pressure under 
which his values for absorbtivity were found. 

In this paper four materials have been chosen to illus- 
trate various relationships. Aluminum is used because of the 
vast amount of previous efforts in the laser weapons program. 
Silver, cadmium and zinc were chosen because of the wide 


variation in their thermal properties. 


22 





τ. THERMAL RESPONSES 


EXC THE MODEL 
After the laser energy strikes the target surface, the 
absorbed energy is converted into thermal energy in the 
material. The thermal energy so deposited in the surface 
is distributed in various ways in the material. The immedi- 
ate result is a sharp rise of the surface temperature. This 
causes a steep temperature gradient in the material which leads 
to heat transport from the surface and into the material via 
heat conduction. Eventually the surface region reaches melting 
temperature and finally vaporization temperature and the sur- 
face melts and vaporizes when enough energy has been absorbed 
to overcome the heat of melting and vaporization. The exact 
development of this 3 phase process and the temperature profile 
into the medium as a function of time is governed by the heat 
conduction equation together with appropriate boundary condi- 
tions. 
The heat conduction equation in three dimensions is: 

oc. 22 = div(K grad T) + A E 

s OE 
where Ps is a solid density, ος is the. heat capacity, K is 
the thermal conductivity and A is the heat produced per unit 
volume per unit time by the laser source. 


4 Tk ¿TA4TKZ/A ας 


A = (1-R)I X T (eq. 3-2) 
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where (1-R) is the absorbtivity, I is the laser power density 
in Watts/cm- and F is the absorbed power density. The term 
4mk/A has the units of enna and is defined as the absorption 
coefficient where k is the imaginary part of the complex index 
of refraction. The value S = 1/4Tk is the skin depth. The 
skin depth is that distance where the intensity of the laser 
is reduced to l/e of its value at the target surface. [In 
metals, this value is on the order of 107 cm for wavelenghts 
in the infrared. 

There are different boundary conditions for each phase 
of the target material. In general even the source term, A, 
is very complicated as absorbtivity varies with temperature 
and phase. The laser pulse is directed at the target at an 
angle and the pulse shape is Gaussian in nature. The target 
thickness effects the temperature profile in the target mater- 
ial. Even the material constants, ος, Oss and K change with 
the wide temperature ranges and phase changes involved between 
room temperature and vaporization temperature. 

In its full generality the problem is quite intractable 
and for the purpose here, significant simplifications must be 
made by adopting a number of assumptions and approximations 
which will now be presented and discussed. 

We will begin the discussion by considering the development 
of the temperature profile for a single solid phase material 
until vaporization temperature is reached at the surface, that 
is, we neglect the complication of the liquid phase. Specifically 


we make these assumptions: 
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l. It willbe assumed that the laser pulse is constant 
with time and the absorbtivity is also constant, possibly 
an effective value calculated by the method of Thomas and 
Musal [14] previously mentioned. 

2. The entire problem will be done in one dimension. 
This one dimensional approach is valid if the laser's beam is 
perfectly uniform over a large focal area. 

3. The materials thermal parameters, Par Cor Ke cbe. 
are assumed to be temperature independent. 

4. Starting temperature for the process will be assumed 
to be at T = 0° so that all temperatures will be given in °C. 
5. The model considered is a laser beam directed at a 
semi-infinite solid target material, figure 12. Under these 
conditions an exact solution of the heat conduction equation 
can be obtained, Carslow and Jaeger [17]. It is also assumed 
that the laser wavelength is in the infrared region where the 

skin depth is very small, on the order of gee GI. 


Under these conditions the heat conduction equation 


reduces to 


2 πα Ζ 
+ " l = SOS a ae (eq 53) 
92 


where k = Κ/ρ.ς. and a is the absorption coefficient 41Tk/A. 


The initial and boundary conditions are T = 0 at z - o and 
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FICURE 12. MODEL SHOWING LASER RADIATION DIRECTED AT 
A SEMI-INFINITE SOLID AND THE CREATION OF 
THE VAPOR FRONT IN DISTANCE Az 
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The solution to the problem is: 


ο. E Z 
Ezt) = πι νκε/π e Eu c z/2v«t] (eq. 3-4) 

where 

erfc = l - erf 
and 

DNE 
erf = ES f e τ dx 
Υπ ο 


This model allows us to estimate the time when the surface 


reaches vaporization temperature. 





2F 
= E O = 
T(0,t,,) = T. 2 Kc ze "Kt /m (eq. 3-5) 
er 
T K^ 
K ES 


Braginskii, et al., [18], shows experimental support for 
this equation. Harrach [19] uses the so-called "heat balance 
integral" method to solve the problem of moving boundaries 


and obtains a similar result for t 


τη Ξ (311 - 28 *]/« HKT,,/F, 3 (eq. 3-7) 


The numerical ratio of the two solutions for t. is 0.99. 
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Bhnesdutrycultv with this solution to the heat conduction 
equation is that it does not yet involve phase changes. In 
order to remove mass from the target surface, first the 
vaporization temperature must be reached, i.e., enough energy 
must be deposited into the target surface to raise the tem- 
perature from T. = 0 to T. = Ty’ or greater (To = surface 
temperature). In addition, enough energy must also be de- 
posited to overcome the latent heats of fusion and vaproiza- 
tion. We now need a simple model for this second stage in 
the process. 

Figure 12 illustrates the model of a semi-infinite solid 
material. As vaporization is reached, there will be an 
instantaneous removal of material from the target surface. 
This results in a moving boundary condition. One would now 
have to solve a heat conduction problem with these moving 
boundaries, where the solid interface (again neglecting the 
liquid phase) moves at a recession velocity U. 

If this velocity can be calculated, then one could deter- 
mine the mass removal rate. The solution of this heat con- 
duction problem with its moving boundary in one dimension has 
been given by Harrach [15] in analytical form using again the 
"heat balance integral" method. An estimate of this velocity 
can be obtained from a simple heat balance expression. If 


Az is the distance that the interface moves in time At then, 
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where L is the total latent heat from room temperature to 
vaporization temperature, L = τ᾿ Co is some overall 
heat capacity and Po is an overall density (ος and o, may 


be assumed to be those values for the solid to simplify the 


problem). Then: 


AZ τ a 
E o Us (eq. 3) 
> 
Us = ο τοπ (eq. 3-10) 
S S V 


This velocity could be thought of as a steady state velocity. 
This expression for the steady state recession velocity of 
the evaporating surface is identical to the one Harrach ob- 
tains from his model. However, in Harrach's model, at times 
close to the vaporization time, the surface recession velocity 
is lower. As no closed form expression exists for times near 
the vaporization time, he shows a plot obtained numerically 
to illustrate the buildup of the recession velocity. This 
solution is reproduced in Figure 13. The result of this 
analysis is that the value for the steady state surface 
recession velocity obtained from the general heat balance must 
be reduced in the region where total laser pulse time is on 
the order of surface vaporization time. 

The considerations discussed above do not take into account 
yet that the material immediately on the target surface may 
be at a temperature higher than the vaporization temperature. 


Anisimov [20] suggests that the surface is indeed at a higher 
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temperature caused by an elevated boiling point created by 
the dynamic and high pressure situation set up at the target 
surface. We will now expand our simple model somewhat so 
that it incorporates the idea that the surface may be at a 
higher temperature than the vaporization temperature and 
that the vapor, when it does leave the surface, leaves at 

a velocity of U, - /kgT_/m. This would be the z-direction 


component of the rms thermal velocity where DE Ἐν 7 + y 2 εν” 


y Z 
and Ko is the Boltzman constant. We make the following 


assumptions for this new model: 


B. THE THEORY OF EVAPORATION 
l. Assumptions 

We will now assume that with the laser energy directed 
at a target an amount of material of depth Az is raised above 
the vaporization temperature and then is removed, Figure 12. 

The materials thermal properties, as before, are 
assumed to be constant with temperature. The liquid and solid 
phase densities and heat capacities are considered to be approxi- 
mately equal and the gas phase heat capacity, Cur is equal 
29. 3/5 Ca" These assumptions are based on data obtained by 
experiment, Kelly [21], CRC Handbook [11], and JANAF Thermo- 
chemical Tables [22]. In addition, it is assumed that the 
material in the small increment Az, to be removed, has heat 
capacity = and solid density, og: as it has not had time 
to expand significantly while being heated to a temperature 


EL T. 
S V 
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2. Temperature Profile in the Material 
Under the above conditions and the one dimension 
model of semi-infinite target thickness, one would expect 
the temperature profile to be as before (eq. 2-4): 


ΞΕ 


2 
T(z,t) = -πξε7τ e^ Mt 


5 erfc[z/2vxt]) 


This expression is cumbersome to work with as it contains 

the error function. However, it can be approximated fairly 
closely. Let D, the diffusion length, be the distance where 
the temperature decays to l/e of its surface value, D 
Solving the above equation, the distance D is found to be 

D = 0.97/«t or simple D - “kt. It is readily apparent that 
the temperature profile is of a decaying exponential form. 

If the diffusion length is used as the spatial decay constant, 
-σ/Ὀ 


the equation T= Te 


= can be compared to the actual solution. 


The comparison, Figure 14, illustrates the close approxima- 
tion. In fact, the difference is only a few percent and even 
less if one is interested in a small distance, Az from the 
surface, Figure 15. 
See energy Balance 
An energy balance equation can be written to describe 


the heat conduction problem using the temperature profile 


oo -z/D AZ -z/D 
F t = pbdz +p Cs xj Te dz+p CT Az +o ce, d (T,e - T,)dz 


(eq. 3-11) 
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In this equation, the terms p LAz and oS, CST, Az are the energy 
required to raise the small increment Az to vaporization 


temperature and overcome the latent heats. The term 


pe f (T_e72/D - T )dz 


is the energy required to raise Az to a temperature T. E T 


The second term, 


eo 
DAC f 1a dz 


NZ 


is the energy lost due to heat transport into the material. 
This balance holds in particular for t = t the laser pulse 


length. Using the substitutions 


e ^z/D 


T. = T. (eq. 3-12) 
D = kt, (eq. 3-13) 
Az = D In(T_/T_) (eq. 3-14) 
C. = 3/5 Co (eq. 3-15) 


the above energy balance equation reduces to: 


C 
ΠΠ ΠΤ: m, "ΠΠ 


F yt 
ΟΡ 


(eq. 3-16) 


Bo 





The parameter F¿Yt, may now be used to calculate the surface 
temperature from this energy equation.  Implicit in this 
derivation is that the pulse duration, t - t. is comparatively 
small while ES is large so that the vaporized material has no 
time to escape during the duration of the pulse and so all 
the energy arriving must be absorbed in the same material 
leading to a temperature rise above the vaporization tempera- 
cure. 

Alternatively, keeping Foto constant one could pro- 
duce a laser pulse where t - t, is very long and Fo small. 
In that case the vaporized material escapes from the surface 
before the pulse terminates. The laser energy heats the next 
surface layer to the vaporization temperature and as it is 
reached the surface layer escapes and so on, so that the sur- 
face temperature stays essentially at the vaporization tem- 
perature. The equation Az = D In (T,/T,,) obviously breaks down 


in this large t_ case as it predicts Az = 0. The equation 


P 
for uos in this situation is 


co 


2 το. E 
a oc ο ο ocn zu leg 3) 


The expression can now be used to calculate Az/Yt, as a 


function of Fo tp: namely: 


. 7 - PE T,YK 


a (Sai a 


S S V 


> 
P 
"O 
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4. The Mass Removed 
The mass removed during the total time, Ens can be 
expressed as Μ-ρ ΔΖ where M is the mass expressed in grams 
per laser focal area. In the case of short pulses we can 


substitute Az = Din (T,/T,). 
M = Pp, "Kt, In (T,/T,,) (eq. 3-19) 
we now express M as a function of the parameter ος. 


ΜΕ. = R 3 Vic ΠΠ (eq. 3-20) 


It must be kept in mind that T. is also a function of Fo tp 
In the case of long pulses, when the surface tempera- 
ture stays close to vaporization temperature we use the 


expression for Az/vt, given above. 


E O E is 
ΜΕ — P„AZF = Fo EUST (eq. 3-21) 


5. Axial Expansion of the Vapor Jet 


For the purpose of the intended application it is of 
importance to estimate the distance which the vapor front 
travels during the duration of the laser pulse. If that dis- 
tance is small, on the order of Az or less, then no significant 
expansion has taken place during the energy deposition and we 


can use the mass removal relation of the high surface temperature 
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model where T. > ας In this case there will be no further 
experimental constraint. One would fire the main discharge 
of the high voltage source when the vapor jet has expanded 
axially across the gap which would be long after the end of 
the laser pulse. 

In the case of long pulses it becomes important how 
far the vapor surface has traveled when the laser turns off 
and all mass has been vaporized. That distance traveled 
must be commensurate with the experimental geometry to give a 
feasible operating condition. 

ΤῈ m is the total pulse time and t. is the time to 
vaporization then (En 7 €, = te is the excess time remaining 
after vaporization in which the mass located in z can be 
removed from the target and travel a distance, 1, from the 
target surface. The vaporized particles will travel at a 


velocity approximately equal to the rms value of velocity 


in a one dimensional Maxwellian distribution. 


U, = /k,T /m (eq. 3-22) 


V 


Schriempf [16] has obtained the same result for vapor velocity 
in slightly different form using simple heat balance equations 


and conservation of mass: 


U = Jy/RT/Z (eq. 3-22') 
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where Z is the atomic weight of the target and R is the gas 
constant, R = 8.314 x o erg °K To The time to vapori- 
zation may be calculated when the expression (eq. 3-18) is 
applied to the time t. when the surface reaches vaporization 


Bp-mperature, T: = T. 


V 
ue T ^ 
O 


using again K = Kp C.- This value is fairly consistent with 
the previously stated values (eq. 3-6) obtained by Carslow 
and Jaeger [17] and Harrach [15] using completely different 
methods, and with the experimental results of Braginskii [18]. 


An expression for the excess time, te, can now be written: 
E a E = t_r —— (eq. 3-24) 


The distance which the vapor surface travels in the axial 
direction during the pulse time is then 


α = O, Ep = t.) (eq. 3-25) 


tevi lecons rider a short” pulse when 2 < Az. For 
the criterion for determining whether we are in the long or 
short pulse region, we can safely use the expression for the 
long pulse which does not involve the elevated surface 


temperature. 
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AZ = m | ρα ο) 
S S V 


Combining the above expression for 2 and Az we have a condition 


for short pulse: 





ΚΤ, Ge Fo ts ο τος 
U. (t. - essen (eq. 3-27) 
VD os ES -- p Ps L +C T3) 


which can be rewritten to give the condition, 


, g? T ^ 
Pe see τ p mi 
E = (eq. 3-28) 
ES F yt =p CT ve 


O D ss v 
Fo Ep! os (b E C. T) | 


If this condition is not satisfied we must calculate the 


value of & where 





K? T^ 
Q = ut, ie) = Us to = πο (eg. 3-29) 
O 
or 
= 2 Ker, 
L = ο ορ) - z ] (eq. 3-30) 


6. Density and Pressure ín the Plume 


Once expressions are found for Az in the two Cases, 
long and short pulse, then we can calculate the density and 


pressure of the vapor plume of the surface. Conservation of 
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mass tells us that 
p U = p_U (eq. 3-31) 


where ar is the vapor density as Az expands into the cavity. 
In the case of short pulses Az = D In(T./T,), where D pc 
E. = AE, - Ἐν) and UL, = yk, T_/m. Substituting, one can find 


the expression for Pp 


Da In(T,/T,) 


Py - Um da — E) 
ac y 
(Gomes) 
Ke T Z 
t = x 
V 2 
K F 
ο 
κι Ps ο ως nl A 
a rea PO 
O V NUES - EK Ty ] 





K 


As before, p /F_ is expressed as a function of Eee 


The pressure, P, is estimated from the dynamic pressure, 


2 
p = y ; 
2 
ar 
F F 
O O 
= [F Vt - /k« ln Td 
EI eu —_—— | (eq. 3-34) 
5 2 — 2 K T 
(E vyt) v ] 
O p c 
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In the case of long pulses, the expression for De; 


is different. As before, 


AZ E O D S S V u 
P 
Here, one can use 
= AZ u 1 a 
Us = Kerzen) = A ER Ten (eq. 3-36) 
Á ο κ. ο) 
P ^t. 


Again using conservation of mass, P Uy = ο one can 


» 


substitute and arrive at the equation for ο... 


[F /t - po C T V«] 
S =e (eq. 3-37) 


ο ους τ ΤΕ ες — 
V S V Gip 


The pressure is found in the same manner as for the short 


pulse and is: 


E o C T Yk] 


P SS V 
τος = U (θα. 3-38) 
Fo V o K?T 2 
5 — 
KF yt 
o p 
Limit on T. to Avoid Ionization 


The question now is what value can T. take in the 
elevated surface temperature model. As pointed out previously, 


if the temperature is allowed to become on the order of 2-3 eV, 


42 





then the ionization degree reaches a critical value where 
the vapor plume begins absorbing the laser's radiation. This 
absorption by the plume occurs at plasma densities given by 


the plasma frequency 


E 2 
2 = ne fem. (eq. 3-39) 
where mo is the free electron density in the plume. If 
electron densities are allowed to rise above this cut-off 

value for propagation of electromagnetic waves, then strong 
absorption will occur in the low-temperature plasma plume. 


This condition for cut-off can be rewritten as 


m Τι: 10 7/5? (eq. 3-40) 
where Ar is now the laser wavelength in centimeters and n, 

is the electron density in electrons per cubic centimeters. 

If we assume one free electron per atom and solid density in 
the vapor at the target surface then Saha's equation, [23], 

can be used to predict the maximum allowable surface temper- 
ature to avoid cut-off degrees of ionization 


n 
ΟΝ LG OOUR/ Te 


e = = 
a = 2-4 x 10 τ. "Ὁ 


where 


n, = neutral particle density 
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U; = ionization potential of the target material 


in eV. 
The degree of ionization is given by n/n, where n, = ης Τη 
and can be written as follows: 
x -11600U./T 

DCN C" qM ec I c 

and 
/ 2 

n. -Ç + JC + 4n,C 

ES E DAS E (eq. 3-42) 

n, zn, 


We now have a rough temperature criterion for significant 
plasma production to occur. Namely, arrange for the value 
-- such that T, will always be below the value required 


Bcrocut-off. 


n . ΜΖ κ 103), i 
e τ P 
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IV. RADIAL EXPANSION OF THE VAPOR PLUME 


As the vapor plume propagates along the axis it will 
also expand radially. To assess the feasibility of the 
approach one needs to know whether the radial expansion is 
commensurate with the expected geometry. Therefore we need 
an estimate of the radial expansion as a function of the 
distance from the target surface. 

A simple estimate can be obtained if we assume that the 
speed with which the plume radius expands is equal to the 
sound speed, see Figure 16. 

We consider therefore a slice of the plume which has 
traveled along the axis for a time, t. For the purpose of 
this discussion we neglect any expansion of the slice in the 
axial direction due to possible pressure gradients in that 
direction. As the slice propagates along the axis and 
expands radially the temperature and density must change and 
so the sound speed is a function of the temperature which is 


a function of the z-position on the axis. 


dR(t) 
che 





= am = YC (Y-I) T (eq. 4-1) 


The expansion is assumed to be adiabatic so that 


Fi 
D 

rn 

d 


(eq. 4-2) 


7 
oO 
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dR/dt = a 


TARGET 


R(t) 





GR dite 


FIGURE 16. MODEL ILLUSTRATING THE RADIAL EXPANSION 
OF THE VAPOR 
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where To’ E. and o. are measured at the target surface and 
Ti, Pie and p]; are taken to be at some later time, after 
the vapor has expanded. The quantity y is the ratio of 
specific heats and is 1.66 for monatomic metal vapor. 

For small changes in density and temperature, dp and dt, 


(ot do" T * dT, v/v-1 
T 


a ) = ( (eq. 4-3) 


Or, 


ot do T+dT,1/y-1 
T 


( (eq. 4-4) 


Using binomial expansion for the temperature term yields 
et leq. 4-5) 
y-1 | 


Or simply, 


Se = (es) E (eq. 4-6) 


Since the density decreases with increasing radius (again 
assuming that the thickness of each slice does not change) 


One can write 


νο = p R? : (eq. 4-7) 
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and 


dR, 2 
=) 


dp 
(1 ο (1 R 


= (eq. 4-8) 


Again using binomial expansion for the quadratic in R and 


substituting for dp/p the above expression yields 


dl 
Y-l 


dT Par - B 
Talo) = (eq. 4-9) 


(1 + 
By assuming the second order (dT) (dR) term to be negligible 


compared to the other terms and rewritting gives, 


2dR 


= _ 1 dat q 
B = PET (eq. 4-10) 


Ro is now taken to be the initial plume radius which is 
the laser's focal radius. To is assumed to be the surface 


temperature. Integrating both sides yields 


A ne) Ñ 
(a = ο) (eq. 4-11) 
ο ο 
DES 
Du 21/2 E 
ES G) (eq. 4-12) 
ο ο 
Smbsicituting vT - a//C, 1-1) and a = dR/dt yields 
(E = _ GR/dt ( (eq. 4-13) 
VoM ID 
O E DT, 
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The term ler is simply UL, the initial vapor velocity. 


Therefore, 








(Ry t-y - i dR 
R oe ay Gls 
O M 
and, 
dR E ES 
l=y | 1-y dt 
R R 
O 








f dR = m 
R P E 
O O 
or, 
YU 
Bey ge V z 
(I) = 1 + p (t Ἐν) 
ο η Ο 


But, UL, = (oat) is the axial distance, £, that the vapor 


plume will travel in the cavity. 


EE e yd 
(3 ) = l + R 
O O 
Or, 
R = R {1 + y2/R e 
O O 
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Therefore, 


(eq. 


(eq. 


(eq. 


(eq. 


(eq. 


(eq. 


4-14) 


4-15) 


4-16) 


4-17) 


22090) 


4-19) 





This expression can now be used to calculate the radial 
expansion as a function of the axial distance traveled 2. 
The previously derived equation for temperature can be 


used to calculate T at that distance R: 


= B T.l/vy-1 
(I = (m) (eq. 4-20) 
Q O 


The density and pressure can be obtained using the adiabatic 


expansion relationship: 


ODA u pe by 
O 
mc xem L/Y-l ^ 
= m πο (eq. 4-21) 


It must be noted again that this calculation gives merely 
an estimate for $ and in a more rigorous model would have 


to be replaced by a proper fluid dynamical model. 
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ο ΡΕ 


There are now two cases of the mass removal process. 
Both involve Fo low enough to avoid significant ionization 
so that all the laser's energy is either absorbed at the 
material's surface or is reflected. There will be little 


absorption in the vapor plume. 


CASE I. SMALL PULSE TIME 


This is the situation where the following condition is 


1s met: 2 2 
2 K > 
U [(F_vt_)~ - ] 
O K 
ος Eee epee eee 
F A V 
Oo p pg(L + CST.) 


In this case the surface temperature will be elevated slightly 
above the vaporization temperature and is described by equation 
(3-28). Figure 17 is a plot of the elevated surface tempera- 
ture, T. as a function of noter the controlling parameter 

at 10 um. On this plot, the lowest value of TS for each 
material corresponds to T, and Ft. The maximum value 

Lor T. is the temperature at which there will be cut-off and 
significant absorption of the radiation occurs by the vapor 
plume due to ionization. That maximum To was calculated by 
first finding the maximum allowed free electron density, nor 
given by equation (3-40). Next, calculate the maximum 


allowed degree of ionization, no/n.. Then, use Saha's 
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equation, equation (2-42) to find the maximum allowed To: 
It will be a different value for each material and wave- 
length, and can vary by several thousands of degrees. 
Figure 18 is a plot of MF, vs. Bares for Case I of 
Short laser pulses and Case II of long time pulses from 
equations 3-20 and 3-21. The difference between Case I and 
Case II is slight except at high values of Fo ts. Case I 
gives a lower value of Με, than Case II. In this plot, the 
maximum value for Case I corresponds to a Fo”tp and a T. 
where plume absorption begins. Case II was also plotted 
to that same value of PEU however, the surface tempera- 
ture does not rise above vaporization temperature and could 


be plotted to much higher values of Foren. 


DEE CASE II. LONG PULSE TIMES 


Case II is the situation where the condition on 


2 p 
U [(F t) - ———] 
V Q K 
F > = 
O Have - pC .T vk 
RENI ] 
ο Ῥ P LC T) 


is not met. This is the case of low ES and long E The 
surface temperature remains at essentially To In this 


case abosrption by the plume never occurs. 


C. MASS REMOVAL CALCULATION 
A numerical example is now presented to show the use of 


the various equations and to illustrate that masses on the 
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order of milligrams could, theoretically, be removed without 
plasma production at a materials surface. The material 
chosen for the example is Zn, as it has a relatively high 
Z-number and low latent heats and vaporization temperature. 


The following thermal properties are from Table II: 


o. - 7.13 g/cm? 
C = 0.39 Joules/g°K 

2 
κ = 0.41 cm /sec 
K = 1.13 Watts/cm *K 
[P = 907 °C 
V 
Z = 65.38 

B 4 1072 2 
Let us assume that Da = 10 Watts sec Tem: and 


that the laser's focal area = 107? ene First the condition 


for the two cases must be tested from equation 3-28, Figure 
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ΠΕ Fo is chosen to be less than 6.22 x 10° then Case II 


must be used. If not, Case I must be used. Let us choose 


B = 2 x MO ote ence Then, 


5 


ct 
N 


(10°/F_)* Bess lo 


25 usec. 


Prom Equation 3-21 or Figure 18 for Case II, 


MES TS ο. Watts grams cm 


4 


or, 


OS grams cm ^ 
If the focal area is assumed to be 107? ent then the mass 
removed is 3.79 mg. The axial expansion during the pulse 


t č 15: 
p 


The density and pressure before radial expansion can be found 


using equations 3-37 and 3-38. 


2.29 x 10 ^ g/cm? 


D 
N 


P = 2.65 x Πο ergs/cm” = 3.85 psi. 
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If the anode-cathode separation is assumed to be 4 cm and 
Ro = 0.18 cm, then the radial expansion at z = 4 cm can 


be found from equation 4-19. 


R= «165 cm. 


The pressure, density and temperature after radial expansion 


to R= 1.65 cm and axial expansion to z = 4 cm are calculated 


from equation 4-21. 


a ee 
Een 
Dee 


This shows that a relatively low power CO, laser pulse 
OT x 10° CET of 25 usec duration can produce a 


zinc vapor plume with the desired mass and dimensions. 
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VI. CONCLUSION 


As stated in the Introduction, the purpose of this study 
was to assess the general feasibility of using a laser's 
radiation to remove masses on the order of milligrams from 
a target's surface. This material will fill a discharge 
channel between the anode and cathode of a high voltage 
pulse device which, when activated, will couple electrical 
energy into the material in the discharge channel causing 
the atoms to undergo electronic transitions and then giving 
off intense x-rays. 

The model presented herein is a heat balance model at 
the target surface. No attempt has been made to solve the 
complicated heat conduction equations with its associated 
moving boundaries. In this model, the controlling parameter 
is Bayt, the absorbed power density and the laser pulse time. 
Two separate cases emerge from the model. Case I is the case 
where, keeping ος constant, ES takes on a relatively high 
value while the pulse time is small. Case II considers 
the situation of the same value for E VE but smaller Fo 
and longer pulse time. Using this model, theoretically, 
masses of the order of milligrams can be removed from 
metallic surfaces with relatively low power densities. 

The mass removal model also predicts low temperatures 
at the target surface and thus in the discharge channel. 


As stated in section II this will produce low degrees of 
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ionization. In order to increase the coupling efficiency 

of the electrical energy into the vapor plume, the ioniza- 

tion level must be increased. This could be done by increasing 
the laser power level towards the end of the pulse or by using 
microwave radiation after the plume has propagated into the 
discharge channel and prior to the application of the high 
voltage pulse. 

In addition to mass removal, the controlling parameter 
Bet is used to predict the temperature, density and 
pressure of the vapor plume as it begins to propagate. 

As the plume propagates into the vacuum chamber it undergoes 
axial and radial expansion. Temperature, density and pres- 
sure decrease during this process and can be predicted using 
a model of adiabatic jet expansion presented in section II. 

The mass removal model predicts that a laser can remove 
masses from a target surface large enough to be used in 
the proposed experiment. In addition, the expansion of the 
vapor plume is sufficiently small as not to damage the com- 
ponents located in the vacuum chamber. Using a long laser 
pulse (Case II) the proposed technique shows that the vapor 
in the discharge channel is relatively cold and requires 
pre-ionization prior to application of the high voltage 


pulse. 
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APPENDIX A 


The interactions of materials with electromagnetic 


waves is formulated on the basis of Maxwell's equation: 


where, 
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We will assume that our concern is with non-magnetic, 
electrically neutral materials. In conducting mediums it 
is the contribution from the free electrons that predominate. 
This will allow us to ignore the polarization effects. 


Maxwell's equations now reduce to the following: 


=- > . 9H 

ο; Ho 5t 

A £ 3E _ 
= —— + 

V x H es ^t d 

VE = 0 

V.H = 0 


The general wave equation for E is then found by taking the 
curl of the first equation and the time derivative of the 


second equation. 


Beau Ze 
ὃ τσ τι.  ὃ 9E IT 
Beer TE 
Using the vector identity 
— Ξ5 = a et M 
V x Vx E = V(V-E) - E 
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and noting that V.E = 0, we can exchange the order of 
differentiation in the electric field equation and multiply 


the magnetic field equation by ug to give 


τ 





τη 
JE = 
ιά 9 (yx H 9 E ðJ 
RES "o nd Po) ox Ho*o p - 
ot 
Adding and substituting C = : r 
CUORE 
ο ο 
2 = INDEM. Nm 
S eee oF ct 


Plane harmonic waves solutions for the electric field are 


of the form 


E s Eo exp{ [-27kz/A] + [iw(t- nz/C) ]} 
where 

n = Real part of index of refraction 

k = imaginary part of index of refraction. 


Substituting this into the electric field equation derived 
above along with Ohm's law, J = oE = -Nev, gives us an 


expression for n-ik. 


n-ik = vl- io/e u 
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g is the conductance and is in general a complex quantity, 


namely σ - σι 7 105. Then, 








n-ik = 1 - 


In order to evaluate n and k one can use the Drude-Lorentz 
model. In this model the atom and its outermost electron 
are assumed to follow the equation governing the motion of 


the harmonic oscillator: 





where 


effective mass 


m* 


relaxation time. 


A 
11) 


EE E cIurion of the form e UF is assumed for Y and E- 
and Ohm's law, cE = -Nev is substituted for v, one can show 
that 

Sere 


m* (1+ wt) 


By separation of the real and imaginary parts one can 


write an expression for σι and PE 
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2 L+ ut! 


Q 
N 


where σ. is the static Of ος conductivity. As the frequency 


approaches zero, w > 0, 


Ne? rt 





The real and imaginary parts of the index of refraction, 


n and k can now be found using a complex identity: 
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T Can now be eliminated from the expression for ση» and σ2 


by using the equation for σα” The resulting equations for 


o. and oa are: 
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If normal incidence is assumed, the abosrbtivity of a 
material can be calculated, again using Maxwell's equations 
and the appropriate boundary conditions. Figure 20 illus- 


trates the geometry. 





= | . 9E 3H 
Par ea + 
EN ο. - EL 
IES aT 

As before, 
Ie c "η ΠΙΕΡΟ} 


Substituting this expression for electric field into 
Maxwell's equations yields a relationship between H and 


E: 
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MEDIUM 1, VACUUM MEDIUM 2, METALLIC SOLID 





FIGURE 20. BOUNDARY CONDITIONS FOR E AND B FIELDS 
AT THE TARGET INTERFACE 
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The subscripts 1 and 2 are for the mediums. Medium 1 is 
vacuum and medium 2 is the target material. The subscripts 
i, r, and t denote incident, reflective and transmitted 
waves. Using these boundary conditions and the relationship 


for Es in terms of E. the ratio of E;/E, can be readily 


obtained: 
Bi 
E. = n-1-ik = E 
τ n-l41ik 
where r is the amplitude reflection coefficient. It has 
been assumed that n, = 1 and k, = 0, which is the case in 


1 È 
vacuum; The reflectivity, R is given by |r * r| and there- 


fore is written as: 


Ban 
(n+1) 2 + k? 


The absorbtivity, A, is the desired result of the derivation 


where: 


A = (1 - R) 


The abosrbtivity can be written as: 


án 
A 2 2 
(n+1)7 + k 
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